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Abstract 



We analyze electromagnetic waves propagation in one-dimensional periodic 
media with single or periodic defects. The study is made both from the 
point of view of the modes and of the diffraction problem. We provide an 
explicit dispersion equation for the numerical calculation of the modes, and 
we establish a connection between modes and poles and zeros of the scattering 
matrix. 

Periodic media with defects have been intensively addressed in the field of Quantum 
Mechanics (see [[j]-|3|] an d references therein) and also, since the development of photonic 
crystals, in Electromagnetics 0-0]. On the general subject of photonic crystals, the inter- 
ested reader can find an impressive bibliography on the Internet Pj. From the theoretical 
point of view, the quasi-totality of the studies are involved in the characterization of the 



spectrum of the infinite medium. Nevertheless, for the working physicists, the main problem 



is that of the finite structure. Indeed, one can only imagine experiments, such as the diffrac- 
tion of a plane wave, by finite devices. The main question is then to relate the modes with 
the behavior of the diffracted field. The simplest connection is that of the determination 
of the conduction bands: whenever the frequency belongs to the spectrum of the infinite 
structure, the finite one allows the guidance of waves while in the gap the electromagnetic 
field decreases exponentially. In the present paper, we study, as a model problem, the simple 
case of a one dimensional periodic structure with defects, which may model for instance a 
quantum cavity. Such a device has also been intensively addressed P- jl^ll - Nevertheless, 
these studies involve the infinite structure for which a spectral analysis is given, whereas in 
the present communication we aim at establishing a link between the diffractive properties 
of a finite structure with a defect, and the known spectral properties of an infinite structure 
with a defect. 

Throughout this paper, an orthonormal triaxial cartesian coordinate system (0, x, y, z) 
is used. We consider a periodic structure, described by a bounded real one-periodic function 
e(x), (e(x) = e(x + 1)), representing the relative permittivity with respect to x, whereas 
the permeability is assumed to be (j,q, i.e. that of vacuum. The structure is assumed to be 
invariant in the y and z directions and the harmonic fields (time dependence of exp (—iuit)) 
are invariant along z. That way, the field is described by a function u n (x) exp(iay), where 
a e (— 7r, +7r] is the Bloch frequency (in case of a scattering problem, the frequency a is 
equal to k sm8, where 9 is the angle of incidence of an incident plane wave). When the 
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electric field E is parallel to the z-axis (E\ | case), Un(x) exp(iay) represents the ^-component 
of E and when the magnetic field H is parallel to the z-axis (H\\ case), it represents the 
^-component of H. 

Denoting: (3 2 (x) = k 2 e(x) — a 2 , /3 2 = k 2 . — a 2 , and setting U = (u, q~ l d x u), the propa- 
gation equation takes the form: 

^=(_ ? V?) C/ (1) 

with: q(x) = 1 for E\\ polarization, q(x) = e (x) for H\\ polarization. The monodromy 
matrix of the equation, or Floquet operator, is the 2x2 matrix T k ,a such that: U(x + 
d) = T ka U (x). When considering only n periods of the medium embedded in vacuum 
and illuminated by a plane wave (the device extends over [0, n]), the following boundary 
conditions hold: 

iPou\ x=0 + d x u\ x=0 = 2i(3 ^ 

if3oU\x=n d x u\x=n 

from which the reflection and transmission coefficients can be derived: 

r n = u\ x=0 - 1 

tn ^|a:=7i 

In the infinite medium, the conduction bands are characterized by the condition \tr (T)| < 2, 
we thus define: 

G = {(k,a)/\tr(T k>a )\ > 2} 
B = {(k,a)/\tr(T kta )\ < 2} 

For (k, a) G GUB, we denote (v, w) a basis of eigenvectors of T^ ^, such that det (v, w) = 

1, associated to eigenvalues fi and — ( as a convention < 1 for (k, a) e G). We write in 

li 

the canonical basis of C 2 : v = (i> 1; v 2 ) , w = (wi, w 2 ). 

Let us now introduce a defect in the crystal in the following way: we replace one period 
of the crystal by a layer of width h and permittivity (x) (e^ is simply assumed to be 
real and bounded). In view of the spectral problem, the structure is therefore made of two 
periodic half-spaces connected by an inhomogeneous layer extending over [0, h}. We denote 
by To the monodromy matrix of the defect, omitting here to explicitly write the dependence 
in (k,a). We have the following results: 

Proposition 1( Pj3|,|i"2]j ) : The defect does not modify the conduction bands. 

Proposition 2: The defect modes of the structure correspond to couples (k,a) such that 
det (T w, v) = 

Proof. As we introduced a defect, we may now accept an increasing solution in the left 
semi-crystal and a decreasing solution in the right semi-crystal. The problem is to match 
these two solutions. This is only possible if T switches increasing solutions to decreasing 
ones, that is if Tow 6 Vect(v) QED. 

Suppose that the structure is finite and that the defect is switched between n periods. 
The monodromy matrix is: T^ a T T'^ Q . In basis (v, w), T is written: 



2 



Denoting: 



a c 
bo d 



i w 2 -if3oWi w 2 + i/3 wi , 
X {Xij) 1 -v 2 + ip Q v x -v 2 -ip vi ' ' 



an expression of the coefficients (r n ,t n ) can be easily obtained: 



where: 



Or \ P^ 2H ) + ft. \ ~ 2 ^0 ^ m 



P PO = X2iXnao^ 2 + (%2i c o - Xu & o) X - X2lXll^0 

? ( X ) = -X2lXl2fl0^ 2 + (-X21X22C0 + XllXl2&0) X + XllX22^0 • 



Assume now that there exists some defect mode for a couple (k ,ao) so that |«o| < &o, 
allowing to define an angle of incidence 9 by a = k sm9 Q . 

It is easily seen that the equation of Proposition 2 simply writes: do (ho, «o) = , whence 
we get from 

/, /, x -Xiibo + xliCo + X2iXn/^ 2n a 

r„ (,fto, ^oj 7 kz — 

X11X12O0 - X21X22C0 - X2iXi2/U^ n ao 

-X21X22C0 + X11X12O0 - X2iXi2/i^ n ao 
As n tends to infinity, we have the following limits: 

r n (k , 6 ) — ^0 " * n (k , 6 ) — ^ 

X21X11O0 - X22X21C0 X21X11O0 - X22X21Q) (5) 

r n (fc, O ) — ► — , for k ^ k t n (k, 9 ) — > 0, for k ^ k 

X22 



Remark 1: Obviously from (|]): 



X21 



- 1, and so as a conclusion \r n (k , 9q)\ tends 

X22 

pointwise towards a value that is strictly less than 1, whereas for any point different from 
k , in a small enough neighborhood of k , \r n (k,9 )\ tends to 1 as n — ► 00. This result 
means that the reflected energy admits a sharp minimum near k Q . It is important to note 
that the minimum of \r n (k, 9q) \ is not a priori reached for value k of the wavenumber, but 
the sharpness is all the more important as n is important. Clearly, for a given incidence 9o, 
it is possible to transmit waves of wavenumber belonging to a small interval near k and 
tending to {k } as n tends to infinity 

It is known that for a fixed incidence 9, the scattering matrix admits a meromorphic 
extension to the complex half-space {Im(A;) < 0} [|17]. From @, poles and zeros of r n are 
solutions respectively of equations: 
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/i 4n a + /i c o = do 

X21 



(6) 

4n X2iXu 2n X21 X12, \ _, 

jj, a + /i Co Oo — ao 

X11X22 X22 / 

and the limit of both equations as n tends to infinity is just do (k, k sin 6) = 0, which is 
the equation defining the defect modes. As both equations (||) are analytic perturbations of 
d (k, /csin 9) = we can conclude: 

Proposition 3: One defect mode (ko, /cosin# ) is associated to exactly one zero k™ ofr n 
and one pole k™ of r n and t n . Moreover k™ and k™ tend to k as n tends to infinity. 

Remark 2: This suggests that, for sufficiently large n, there exist two positive numbers 
7„ and 5 n such that: 

k n z =ko- ij n S n + O (S n ) 
k n v =ko- iS n + O (5 n ) 

We know that, as n tends to infinity, S n tends to and that r n (ko, 6o) admits a limit 

given by In the vicinity of ko, we can write r n (k, 6) ~ whence we obtain: 

X22 k - k™ 

n a x X2iIm(A;^) X?iX 2 i 1& o - X21C0 

r n {ko, ) ~ — v— y, so that 7„ 



X22 Im {kfj ' " X21C0 - X2iXnX 2 2 1& o ' 
Remark 3: For real z, z 1— > is the equation of a circle of diameter \/l + y%, 

X22 Z — K p 

so that while in the gap the reflection coefficient belongs to the unit circle of the complex 
plane, it describes a circle for real k varying in the vicinity of /co- 
in order to be able to apply Bloch- waves method to defect structures, an approximate way 
consists in periodizing the defect, it is the so-called supercell approximation. Let us apply 
this method to our simple example. We consider thus an infinite periodic structure whose 
period is made of one defect switched between n periods of the previous medium. We call 
it a super-structure. We know that the global monodromy matrix of the super-structure is: 

tr (t) < 2. Now 



Tk, a = T^ToT'^, and that the conduction bands are characterized by 
what happens at k = k ? Obviously:tr ^Tk,a j = f 211 ^ + ^~ 2n do so that for (k, a) = (k , ao) 
: tr (^Tk ,a j = A i2 " a o an d thus tends to as n tends to infinity. This means that, at a = ao, 



there exists an interval J n of wavenumbers k over which tr yTk,a J < 2 and therefore J n is 
a conduction band in the super-structure. Of course the length of this interval depends upon 
n and tends to as n tends to infinity because J n tends to {k }. This is a crucial point as 
it justifies the use of Bloch-waves theory in the framework of the super-cell approximation, 
at least in a small neighborhood of ko- Finally, we can state: 

Proposition 4: In the super- structure, the existence of a defect (ko, ao) implies the 
opening of conduction bands inside gaps of the unperturbed structure that can support defect 
modes. The widths of these " defect" -conduction bands tend to zero exponentially with the 
number n of subperiods constituting one period of the super- structure. 

Remark 4: For a finite superstructure, and as the medium is periodic, we may apply 
a theory [|TTJ allowing to compute the superior envelop of the modulus of the reflection 
coefficient of a periodic structure. In the present problem, this is just the graph of: 
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where T 
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We have analyzed wave propagation in one-dimensional photonic crystals with one defect. 
We have shown a connection between the scattering properties and the modes: defect modes 
of the infinite structure give rise to a pole and a zero explaining the behavior of the refection 
coefficient. These results should help in the understanding of more complicated structures 
such as bidimensional photonic crystals, for which our present work could be used as an 
approximate theory HQ- 
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